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We investigate the scattering phenomena in two dimensions produced by a general finite-range 
nonseparable potential. This situation can appear either in a Cartesian geometry or in a het- 
erostructure with cylindrical symmetry. Increasing the dimensionality of the scattering problem 
new processes as the scattering between conducting channels and the scattering from conducting to 
evanescent channels are allowed. For certain values of the energy called resonance energy the trans- 
mission through the scattering region changes dramatically in comparison with an one-dimensional 
problem. If the potential has an attractive character even the evanescent channels can be seen as 
dips of the total transmission. The multi-channel current scattering matrix is determined using 
its representation in terms of the R-matrix. The resonant transmission peaks are characterized 
quantitatively through the poles of the current scattering matrix. Detailed maps of the localization 
probability density sustain the physical interpretation of the resonances. Our formalism is applied 
to a quantum dot in a two-dimensional electron gas and to a conical quantum dot dot embedded 
inside a cylindrical nanowire. 
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I. INTRODUCTION 

There is a permanent requirement of shrinking the semiconductor devices in integrated circuits 1 -. As feature sizes 
shrink into the nanometer scale regime, the device behavior becomes increasingly complicated since new physical 
phenomena at short distances occur and limitations in material properties are reached. In order to keep the good 
characteristics for transistors, new transistor architectures were developed progressively in the last decade. 

Nowadays, there are developments of planar MOSFETs (metal oxide semiconductor field effect transistors}^ as well 
as of gate-all-around (GAA) MOSFETs^^. Both systems are also strongly related to more fundamental research 
structures developed in the last years, like in-plane-gate transistors^, single-electron transistors^, silicon-on-insulator 
planar double-gate transistors^, non-planar double-gate FinFETs^, non-planar trigate transistors 1 ^, nanowire-based 
field-effect transistors (FET)^ 1 -, nanowire resonant tunneling diode o 12 ' 13 , nanowire lasers 1 ^, or nanowire qubits 1 ^, 
whose maturity has still to be proven for industrial applications. Their structural complexity has also progressively 
increased, allowing for double-barrier structure a 12 ' 13 , or multiple core-shell layer a 14 d 6 . The material composition 
includes mainly III-V materials GaAs/AlGasAsSi 1 ^, InAs/InRi 2 ., GaN/InGaN^, but also group IV materials Si 4 ^, 
and Si/G e 11 ! 15 , predominant in the industry. 

The transport phenomena in these mesoscopic devices go beyond the semi-classical limit, and a quantum mechanical 
description of the current and charge densities 1 ^ is necessary. The most appropriate method for analyzing semicon- 
ductor devices with an active region in the nanometer scale and which are almost open (i.e. showing a strong coupling 
between the active region and contacts) is the scattering theory. 

This paper is focused on systems for which the scattering process is a two-dimensional one. Such systems can appear 
either in a Cartesian geometry, like for devices tailored in a two-dimensional electron gas (2DEG)&£, or in a cylindrical 
geometry, like for nanowire-based de vices 3 -^ ^ 1 2 1 1 3 . In these systems there is a strong confinement of the motion in one 
direction, called transversal direction, while the transport occurs in the other direction, called longitudinal direction. 
The scattering problem is a two-dimensional one because the scattering potential is nonseparable, and also the 
incoming electrons can choose different energy channels for transport which are mixed due to the scattering. 

We present in this work a general method, valid within the effective mass approximation, for solving the two- 
dimensional (2D) Schrodingcr equation with scattering boundary conditions. Its solutions are found using the scat- 
tering theory and the R-matrix formalis m 19 i 20 i 21 ' 22 ' 23 i 24 i 25 ' 26 i 27 ' 28 i 29 . This method is a semi-analytical one, and it 
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gives the scattering functions in each point inside and outside the scattering area and for each energy as a function 
of the solutions of the Wigner-Eisenbud problem. As known from the nuclear physics, the Wigner-Eiscnbud prob- 
lem is the eigenvalue problem of the Hamilton operator for the closed counterpart of the considered open quantum 
system 3 ^. The R-matrix formalism is not only numerically very efficient, but it is also suitable for higher dimensional 
nanostructures with complex geometr y 20 ' 21 ' 22 ' 2 ' 25 and general nonseparable scattering potentia l 31 ' 32 . It can also deal 
with more than two terminal o 26 ' 29 . 

Using the scattering functions we analyze further the transport properties of the open quantum structures, especially 
the conductance. Besides the low dimensionality of such systems, the open character is also an essential feature which 
controls the transport phenomena through the structure. When the quantum system becomes open, its eigenstatcs 
yield resonance states which do not have an infinite life time anymore and which are not strictly localized inside the 
quantum system. In this paper we identify the signature of the resonance in each conductance peak and study the 
influence of the nonseparable character of the potential on the resonances and on the conductance through the system. 

An interesting effect in a multi-channel scattering problem is that as soon as the potential is not separable anymore, 
the channels get mixed. If furthermore the scattering potential is attractive, then it leads to unusual scattering 
properties, like resonant dips in the transmission coefficient just below the next channel minimum energy. As it was 
shown analytically for a S scattering potential^ 3 , and later on for a finite-range scattering potentia l 34 ' 35 the dips are 
due to the quasi-bound-states splitting off from a higher evanescent channel. So that evanescent channels can not 
be neglected when analyzing scattering in two- or three-dimensional quantum systems. These findings were recently 
confirmed numerically for a Gaussian-type scatterer— and also for a quantum dot or a quantum ring 3 -^ embedded inside 
nanowires tailored in a two-dimensional electron gas (2DEG), or inside cylindrical nanowires 3 - 1 -. The high resolution 
maps for the "near field" scattering wave functions presented in Refsi 31 ' 36 ' 37 show explicitly increased localization 
probability around the scatterer for energies of the quasi-bound states, in agreement with the resonant reflection or 
resonant back-scattering interpretation of these dip a 33 ' 34 . The Cartesian and cylindrical geometries present different 
"selection rules" for the intcrsubband transmission^. 



II. MODEL 

The electronic states in mesoscopic systems are easily described within the effective mass approximation whose 
validity requires that the envelope function ^>{E,r) must be slowly varying over dimensions comparable to the unit 
cell of the crystal 3 ^. 

In the spherical effective mass approximation, the envelope function associated to the energy E 3 ]j satisfies a 
Schrodinger-type equation 

(f) = E 3D V{r). (1) 

The so-called scattering potential V(f) contains the information about the confinement in the transversal direction, 
and inside the allowed area it is a sum of the heterojunction conduction band discontinuities, the electrostatic potential 
due to the ionized donors and acceptors, the self-consistent Hartree and exchange potentials due to free carriers, and 
external potentials. We use the symbol m* to denote the effective mass of the electrons, while m denotes the magnetic 
quantum number. 

For systems tailored in the 2DEG, the growth direction is chosen the z-direction, while the plane of the 2DEG 
is (x, y). The wave function in the z-dircction, £(z), is taken as known (the simplest form is provided by Fang and 
Howard; 3 - 9 -) so that the three-dimensional (3D) wave function can be written as 

Hf(E 3D ;r)=£(z)i;(x,y). (2) 

The total energy 

E$D = E 2 DEG + E, (3) 

where Eibeg is the energy of the 2DEG level, and E is the energy associated with the motion in the plane of the 
2DEG. The 3D Schrodinger type equation reduces to a two-dimensional Schrodinger equation 3 - 2 . 

ip(x,y) = E ip(x,y), x e (-co, oo), y € [-d y , dy]. (4) 

In the transversal direction the electron motion is limited at the interval [— d y ,d y ] by a confining potential that we 
have considered as infinite. The nonseparable potential V(x, y) varies strongly with the position only inside a small 
domain (|x| < d x , \y\ < d y ) which is usually called scattering region and is quasi- const ant outside this domain. 
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For cylindrical nanowircs, the azimuthal symmetry suggests to use cylindrical coordinates, with z axis along the 
nanowire^i. As long as there are not split gates on the surface of the nanowire, the potential energy V(r) is rotational 
invariant 

V(r) = V(r,z) (5) 

and nonscparable inside the scattering region. A scattering potential which does not explicitly depend on the azimuthal 
angle 9 imposes the eigenfunctions of the orbital angular momentum operator L z as solutions of Eq. |T]) 

y m (E aD ;r, 6, z) = CmWmfa z), (6) 

where 

p irn9 

CM = P) 



and m = 0, ±1, ±2, ... is the magnetic quantum number. This is an integer number due to the requirement that the 
function e %mB should be single- valued. The functions ip m (r,z) are determined from the equation 



d Id m 2 d 
2m* \ dr 2 r dr r 2 dz 



+ V{r,z) 



ip m (r,z) = Eip m (r,z), r G [0,R], z G (-oo, oo), (8) 



where E denotes here the kinetic energy associated with the 3D motion of the electron inside the nanowire, E = E%r). 
We have also considered an infinite potential outside the nanowire. In such a way, every magnetic quantum number 
m defines a two-dimensional (2D) scattering problem. Furthermore, these 2D scattering problems can be solved 
separately if the scattering potential is rotational invariant. How many of these problems have to be solved, depends 
on the specific physical quantity which has to be computed. 



A. Scattering problem for two dimensions 

We consider the following Schrodinger type equation in two dimensions, (x\\,xj_), denoting generically the longitu- 
dinal and the transversal direction, respectively, 



2 



2m 



■&x h x± + V(x||,a;j.) 



ijj(E;x\\,x±) = Etp(E;x\\,x±), x± G Q,x\\ € (—00,00). (9) 



One could consider here different effective masses in the longitudinal and transversal directions and also, for a layered 
heterostructure, a position-dependent effective mass. These effects can be incorporated in the formalism, but within 
this paper we neglect them for the simplicity of the exposure. 

For the Cartesian geometry^, in comparison with Eq. @ we have 



x,x± = y,Q = [-d y ,d y ], 
d 2 8? 
dx 2 dy 2 ' 
E = E^u — E 2 deg, 



A *v^=— + —> ( 10 ) 



while for the cylindrical geometryai, in comparison with Eq. (|8J) we have 

xh = z, x± — r, fl = [0, R], 



d 2 1 d m 2 d 2 



E = E- 



3D, 



which depends on the magnetic quantum number m. 

The analogy between the both geometries appears more evident considering for the cylindrical geometry the unitary 
transformation U : L 2 ([0,R] x R,rdrdz) — > L 2 ([0,R] x M, drdz), with Uf(r,z) = g(r,z) = ^/rf{r,z). The inverse 
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transformation is W : L 2 ([0, R] x R, drdz) -»■ L 2 ([0, R] x R, rdrdz), with J7 t .g(r, z) = f(r, z) = (l/y/r)g(r, z). In such 
a way, the Schrodingcr operator becomes 



UH^zU 1 = U 
h 2 



-A z , r + V(z,r) 



2m* 



dr 2 



2m" 

to 2 -1/4 <9 2 
r 2 <9z 2 



+ V(z,r). 



(12) 



In turn, the term which contains 1/r 2 and which is specific for the cylindrical geometry plays the role of a potential, 
and the Laplace operator appears as known for two dimensions in Cartesian coordinates. The current scattering 
matrix and the localization probability distribution density of an electron in a scattering state remain unchanged 
under this unitary transformation. 
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FIG. 1: The generic geometry of the 2D scattering problem. 



Due to the localized character of the scattering potential it is appropriate to solve Eq. © within the scattering 
theory The potential energy which appears in Eq. ^ has generally two components: 



V(x\\,Xx) = Vjl(xjl) + V acat t(x\\,Xx). 



(13) 



The first one, Vx(xx), describes the lateral confinement of the electrons and is translation invariant along the parallel 
direction xu. We consider a hard wall potential 



V x (x ± ) 



0, x ± eft 

00, Xx & 



(14) 



which defines a quantum wire for the Cartesian geometry and a cylindrical nanowire for the second considered 
geometry. A parabolic wall like in Ref i 36 ^ 37 may also be considered. 

The scattering potential energy inside the nanowire, V sca u{xn, Xx), has generally a nonscparablc character in a 
domain of finite-range and is constant outside this domain. We consider here the nonseparable potential localized 
within the area f2 x [— dn, dn], see Fig. [U 



G il,x\\ < —d 



Vi, 

V sca tt(x\\,Xx) = { W(x\l,Xx), Xx<Ett 

V~2, xx G 0, x\\ > d|| 



II 

du < xii < dii 



(15) 



There are not material definitions for the interfaces xii = ±dii . Usually, they are chosen inside the highly doped regions 
of the heterostructure characterized by a slowly variation of the potential in the longitudinal direction, practically by 
a constant potential. These regions play the role of the source and drain contacts. 



B. Scattering states 

In the asymptotic regions, > du i.e. source and drain contacts, the potential energy is separable in the 
transversal (i.e. confinement) and the longitudinal (i.e. transport) direction, i.e. V{x\\,xx) = Vx{xx) + V s , s = 1,2, 
and Eq. ([9]) can be directly solved using the separation of variables method 

il)(E;x\\,Xx) = </>(xx)<p(x\\). (16) 
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The function <f>(x±) satisfies the transversal equation 

h 2 



2m 



■Ax ± + V±(xx) 



<j>(xx) = E ± 4>(x_l), 



X\ £ Q 



where 



A*. 




dr 



Cartesian geometry 
cylindrical geometry 



(17) 



(18) 



The hard wall confinement potential requires Dirichlct boundary condition at the boundaries dfl of the interval O, 
4>(dQ) =0. As a remark, for the cylindrical geometry, the boundary r = is an artificial one introduced in order to 
use the cylindrical symmetry. At this boundary it is sufficient that 4>(x±) remains finite. 

Due to the electron confinement in the transversal direction x± the solutions of Eq. (fTT)) define the transversal 
modes, (f> n (x±), with the corresponding transversal energies E± n , n > 1. The cigenfunctions 4> n (x±) depend on the 
geometry (Cartesian or cylindrical) and on the confinement potential. In the case of a hard wall confinement, the 
transversal modes are given for the Cartesian geometry by sine functions^, while for the cylindrical geometry they 
are expressed in terms of the Bessel functions of the first kind^I. The transversal modes form an orthonormal and 
complete system of functions. 

The function tp(x\\) satisfies the one-dimensional Schrodinger type equation called longitudinal equation 



2m* dx'i 



ip(x\\) = (E - E±)<p(x\\), x\i £ (-oo, ~d\\) U (d||,oo), 



(19) 



where s = 1 stays for the source contact (xii < —du) and s = 2 for the drain contact (atii > dii). In the case of different 
effective masses in transversal and longitudinal direction, one can use the corresponding effective mass in each of the 
above equations. 

Every transversal mode together with the associated motion on the transport direction defines a scattering channel 
on each side of the scattering area. The scattering channels are indexed by (sn), n > 1, s = 1,2 for each E. In 
contradistinction to the Cartesian geometry, in the case of a cylindrical geometry there is a set of 2D scattering 
problems, indexed by the magnetic quantum number m, that have to be solved. Consequently the scattering channels 
should be also indexed by m. For simplicity we omit the index m in this section, but we keep in mind that we solve 
here a 2D scattering problem and obtain the scattering functions for a fix value of m. 

If the total energy E and the lateral eigencnergy E± n are fixed, there are at most two linearly independent solutions 
of Eq. (|19|) . In the asymptotic region they are given as a linear combination of exponential functions 



fsn(x\\) 




D — }fc 2 „X|| 



x \\ < - d \\ 
x \\ > d \\ 



(20) 



where A s , B Sl C s and D s are complex coefficients depending on n and E for each value of s = 1, 2. The wave vector 
is defined for each scattering channel (sn) as 



(21) 
(22) 



k sn (E) = koy/{E-E ±n -V a )/uo, 
where ko = 7r/2dii and uq — f^k^jlm* . In the case of the conducting or open channels 

E - E ±n -V s >0, 

k sn are positive real numbers and correspond to propagating plane-waves. For the evanescent or closed channels 

E - E ±n - V s < 0, (23) 

k sn are given from the first branch of the complex square root function, k sn = i\k sn \, and describe exponentially 
decaying functions away from the scattering region. Thus, the number of the conducting channels, N S (E), s = 1,2, 
is a function of energy, and for a fixed energy E this is the largest value of n, which satisfies the inequality (|22[) for 
given values of s. 

Each conducting channel corresponds to one degree of freedom for the electron motion through the system and, 
consequently, there exists only one independent solution of Eq. @ for a fixed channel (sn) associated with the 
energy E, ipf?' (E; xii, x±). For describing further the transport phenomena in the frame of the scattering theory it 
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is convenient to consider this solution as a scattering state, i.e. as a sum of an incoming component on the channel 
(sn) and a linear combination of outgoing components on each scattering channel. One can write the scattering wave 
functions in a compact form 3 - 2 - 

S sie ^(^+ d iO ( j )n ( x± ) + £ S M , sn (E)e- ik ^'^+ d lOcf> n ,(xx), x {] < -d„ 

n'=l 

The step function 9 in the above expressions, with 9(x > 0) = 1 and 0(x < 0) = 0, assures that the scattering 
functions are defined only for the conducting channels. Writing explicitly the position of the interfaces ±<i|| at the 
exponent has advantages for the analytical treatment of the scattering proble m 22 ' 40 . As it is discussed in Refsi 26 ' 31 , 
it is necessary to consider the sum until infinity in the second term of the above expression, in order to keep the 
mathematical completeness of the transversal channels. 

The physical interpretation of the expressions (|24|) is that, due to the nonseparable character of the scattering 
potential, a plane- wave incident onto the scattering domain is reflected on every channel - open or closed for transport 
- on the same side of the system and transmitted on every channel - open or closed for transport - on the other 
side. The reflection and transmission amplitudes arc described by the complex coefficients S sn ^ sn and S s i n i tSn with 
s 7^ s', respectively, and all of them should be nonzero. These coefficients define a matrix with Ni(E) + N2(E) infinite 
columns. For an elegant solution of the scattering problem we extend S(E) to an infinite square matrix and set at zero 
the matrix elements without physical meaning, S s ' n i >S n{E) = 0, n > N S (E), s = 1,2. In this way we define the wave 
transmission matrix or wave-function amplitudes matrioj^. It is also called generalized scattering matrix^. This is 
not the well-known scattering matrix (current transmission matrix) whose unitarity reflects the current conservation. 
The generalized scattering matrix is a non-unitary matrix, which has the advantage that it allows for a description of 
the scattering processes not only in the asymptotic region but also inside the scattering area. 

The three-dimensional scattering states, solutions of Eq. ([1]) can be now written as 

¥ n s \E;x h x ± ,x 3 ) =io(x 3 )^ n s HE;x h x^, (25) 

where w(xa) stays for £(z) in the case of the Cartesian geometry and for Cm(#) in the case of the cylindrical geometry. 
Being eigenfunctions of an open system, the scattering states are ortho-normalized in the general sensed 

dfl [ dx\\ ipl l s) {E;x\\,x±)ifj^,\E';x\\,x±)* = S ss >6nn> ^ (26) 

J-oc 9sn(&) 

where g sn (E) = m* /[h 2 k sn (E)] is the ID density of states. We have to mention that for the Cartesian coordinates 3 ^, 
the measures are dfi = dy, dx\\ = dx, while for the cylindrical geometry- 3 - 1 - they are dQ — rdr, dx\\ = dz. 



C. R- matrix formalism for two dimensions 



The scattering functions inside the scattering region are determined using the R-matrix formalism, i.e. 
they are expressed in terms of the eigenfunctions corresponding to the closed counterpart of the scattering 
proble m 19 ' 20 i 21 ' 22 i 23 ' 24 i 25 ' 26 ' 29 . In our opinion this is a more appropriate method than the common mode space ap- 
proach which implies the expansion of the scattering functions inside the scattering area in the basis of the transversal 
modes <fi n (x±). As it is shown in Ref . 35 , 42 the mode space approach has limitations for structures with abrupt changes 
in the potential or sudden spatial variations in the widths of the wire; it breaks even down for coupling operators 
that are not scalar potentials, like in the case of an external magnetic field. In the R-matrix formalism the used basis 
contains all the information about the scattering potential, and this type of difficulties can not appear. 

Thus, the scattering functions inside the scattering region are given as 

oo 

ii s \E; Xll ,x ± ) = J2 a ln( E )Xl(x\\,x ± ), (27) 
i=i 

with x± £ n and X|| G [— 

The so-called Wigner-Eisenbud functions, Xi( x \\ > x ±)i firstly used in the nuclear physic o 30 ' 43 , satisfy the same 
equation as ipn \ x \\ > x±), Eq. (|9|), but with different boundary conditions in the transport direction. Since the 
scattering function ipn\x\\, x±) satisfies energy dependent boundary conditions derived from Eq. (|24[) due to the 
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continuity of the scattering function and its derivative at x\\ = ±<i||, the Wigner-Eisenbud function xi{ x \\ ? x .l) nas to 
satisfy Neumann boundary conditions at the interfaces between the scattering region and leads 



dxi 
dx\\ 



0, 1>1. 



(28) 



=±d,i 



The hard wall confinement potential requires Dirichlet boundary condition at dQ also for the Wigner-Eisenbud 
functions. xi{d^l,x\\) = 0. As already mentioned for the scattering states, for the cylindrical geometry it is sufficient 
that the Wigner-Eisenbud function remains finite at r = 0. The functions xu I ~> E build a basis which verifies the 
orthogonality relation 



and the closure relation 



/ 



r/O 



/ dx\i Xl(xx,x\\)xi'(xx,z\\) 

J-dn 



(29) 



(30) 



Note that for the cylindrical geometry S(x±_ — x'j_) in the relation (|30|) means 5(r — r )/r. The corresponding eigenen- 
crgics to xi are denoted by Ei and are called Wigner-Eisenbud energies. Since the Wigner-Eisenbud problem is 
defined on a closed volume with self-adjoint boundary conditions, the eigcnfunctions xi an d the eigcncncrgies Ei can 
be chosen as real quantities. The Wigner-Eisenbud problem is, thus, the closed counterpart of the scattering problem. 

In the case of the one-dimensional system without spherical symmetry, it was recently proven mathematically 
rigorous that the R-matrix formalism allows for a proper expansion of the scattering matrix on the real energy axis2^. 
In this section we present an extension of the R-matrix formalism for 2D scattering problem. 

To calculate the expansion coefficients af^(E) we multiply Eq. @ by Xl{ x \\ > x ±) an d the equation satisfied by 

the Wigner-Eisenbud functions by i})n\E\ x\\ , x±). The difference between the resulting equations is integrated over 

fl x [— with the corresponding measures, and one obtains on the right-hand side the coefficient afJ{E). After 

using the Green's theorem and the boundary conditions one finds a\^{E) and feeds in it into Eq. (|2"T]) . So, the 
scattering functions inside the scattering region (x\\ <G [— dikdn], x± £ fl) are obtained in terms of their derivatives at 
the edges of this domain, 



1>P(E;x h x ± ) = ±- [ dx' ± 
ko Jn 



d^iE-x'^) 



R(E;-d\\,x'j_,x\\,x_i) 
- R{E; dii , X_l,x\\ , x±) 



dx' H 



d^ n s) {E;x'x' ± ) 



dxi 



(31) 



where the i?-function is defined as 



Df „ , , x tio^Xj(a:||,a:±)xi(^|,a:±) 
ii(£;;a;||,a;j.,a;||,a;j_ ) ) = — 2^ E - E t ' 



(32) 



i=i 



The functions dipn /dx\\ at xm = dbdn are calculated from the asymptotic form (|24[) based on the continuity conditions 
for the derivatives of the scattering functions on the interfaces between the scattering region and leads. 

With these results the scattering functions inside the scattering domain are expressed in terms of the wave trans- 
mission matrix S 



^{E;x\\,xx) 



/2tt 



0(E)[1 - S T (E)}K(E)R(E;x\\,x ± ) 



(33) 



where the component (sn) of the vector ^> is the scattering function tpn (E; xu , x±), n > 1, s = 1,2 and S T denotes 
the matrix transpose. The diagonal matrix K has on its diagonal the wave vectors (|2ip of each scattering channel 



fT? s k sn (E) 

'{n) = — ; o nn >o ss i, 

k 



(34) 



8 



i' > 1, s,s' = 1,2, and the vector R(E\ x\\ , x±) is given as 



1=1 



where \i nas the components 



(Xi )sn 



Xi(aJ±,(-l) s rf||)^n(a;-L)^, 



(35) 



(36) 



n > 1, s = 1, 2. The diagonal (9-matrix, s n,s'n'(E) = 9(N S (E) — n) 5 ss i S nn i, n > 1, s = 1, 2, assures non-zero values 
only for the scattering functions corresponding to the conducting channels. 

Using further the continuity of the scattering functions on the surface of the scattering area and expanding 
R(E; ±<i|| , x±) in the basis {4>n(x±)} n>1 we find the relation between the matrixes S and R 



S(E)= l-2(l + iR(S)K( J B)) _1 0(E) 



(37) 



with the i?-matrix given by means of a dyadic product 



R(E) = u a J2 



1=1 



Xi Xt 
E-Ei 



(38) 



According to the above relation, R is an infinite-dimensional symmetrical real matrix and its elements defined by 
Eq. (|3"5)) are dimensionless. The above form allows for a very efficient numerical implementation for computing the 
i?-matrix. 

The expression (|37[) of the S-matrix in terms of the i?-matrix is the key relation for solving 2D scattering problems 
using only the eigenfunctions and the eigenenergies of the closed quantum system. They contain the full information 
about the scattering potential and carry it over to the i?-matrix. The matrix K characterizes the contacts and can 
be constructed using only the information about the potential in these regions. On the base of Eq. (|37[) the wave 
transmission matrix is calculated and after that the scattering functions in each point of the system are obtained 
using Eqs. (j2"4")l and (|3"3"|) . Further on, each transport property of the open quantum system can be derived from the 
scattering function in terms of the scattering matrix. 



D. Reflection and transmission coefficients 



Using the density current operator 

J(r) = ^ (*(r)W(rT - *(f)*V*(f)) , (39) 

one can define, as usually, the transmission and reflection probabilities^. Here *S?(r)* denotes the complex conjugate 
of the scattering wave function (|25|) . 

The transversal component of the density current jj_(xj_, xt\, X3) is zero in leads, because 4> n (x±) are real functions. 
The component X3 of the incident density current is also zero, either due to the confinement in the third direction, 
like in Cartesian geometry^, or due to the symmetry reasons like for the cylindrical geometry^. What remains is 
the longitudinal component of the particle density current ju (xj_, xti , X3), which provides after the integration over 
the cross section of the lead with the corresponding measure, d£l, the very well-known relations for the transmission 
and reflection probabilities. The probability for an electron incident from the source, s = 1, on the channel n to be 
reflected back into the source on the channel n' is 

R'nl^J^lSLln'l 2 , (40) 

and the probability to be transmitted into the drain, s = 2, on the channel n' is 

^U^l|5f n , 2n ,| 2 . (41) 
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The reflection and transmission probabilities for the evanescent (closed) channels are zero. The total transmission 
and reflection coefficients for an electron incident from reservoir s = 1 are defined as 

^ (1) =E^" K {1) =Y. R nl- (42) 

More detailed properties of the many-channel tunneling and reflection probabilities are given in Ref.— , but note that 
our indexes are interchanged with respect to the definitions used there. 



E. Current scattering matrix 

Further, we define the energy dependent current scattering matrix as 

S(E) = K 1 ^(E)&{E)S(E)K- 1 ^(E), (43) 
so that its elements give directly the reflection and transmission probabilities 

ln (£)| 2 = rW(E), \S 2n ,, 2n (E)\ 2 = R { nl(E), 
|S 2 „m„(£)| 2 = T<£,(E), \S W , 2 „(E)\ 2 = T3(E)- 

The diagonal 0-matrix assures that the matrix elements of S arc nonzero only for the conducting channels, for 
which the transmitted flux is nonzero. Using the i?-matrix representation of S, Eq. (|37|) . we find from the above 
relation 

S(E) = 0(E) [1 - 2(1 + if2(E))- 1 ] 0(E), (45) 
with the infinite dimensional matrix fl 

OO 

f2(E) = K 1 ^(E)R(E)K 1 / 2 (E) = Uo £ J^L_ (46) 

i=i 1 

and the column vector 

d l (E) = K 1 / 2 (E)xi, (47) 

with I > 1. 

Further we express the total tunneling coefficient in terms of the current transmission matrix, 

T(E) = Tr[a(E)^(E)}, (48) 



where er denotes the part of S which contains the transmission, cr nn > (E) = S-2n' ,in(E) , n = 1, Ni(E) and n 1 = 1, N2(E). 

According to the definition (|46|) the matrix fl is a symmetrical one, fl = fl T , and from Eq. (|45|) it follows that 
S also has this property, S = S T . On this basis one can demonstrate that the tunneling coefficient characterizes one 



pair of open channels irrespective of the origin of the incident flux T^J, 



S-2 



Sir: 



= T$l This is a 

well-known property of the transmission through a scattering system and it shows that the current scattering matrix 
used here is properly defined. The restriction of S'-matrix to the open channels is the well known current scattering 
matri x 20 ! 22 ' 23 , commonly used in the Landaucr-Buttiker formalism. For a given energy E this is a (N1+N2) X (N1+N2) 
matrix which has to satisfy the unitarity condition, according to the flux conservation. 

In the numerical computations, the matrixes S, R, fl, S and have the dimension 2N x 27V, and the vectors \i> 
ai(E) have 2N components, where N is the number of scattering channels (open and closed) taken numerically into 
account. The number of the Wigner-Eiscnbud functions and energies computed numerically establishes the maximum 
value for the index I. 



Resonances 



The relation (|4"5"]) is the starting point for a resonance theory of the transmission through a structure with a 
scattering regio n 22 ' 32 . The singularities of the current scattering matrix S which satisfy the equation 

det[l + iS2(E)}=0 (49) 
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are usually classified as bound states and resonances. The bound states are characterized by real negative energies 
while the resonance energies, E i = E i — &T;/2, I > 1, lie in the complex energy plane below the real positive axis 
according to the causality^. The scattering matrix S and consequently the total transmission T(E) are defined only 
for energies in the continuum spectrum (E real positive) of the scattering problem and they are analytical functions 
over the whole domain. Although they have no singularities in the definition domain, their energy dependence is 
determined by the resonances, especially by those ones which lie in the vicinity of the real axis. In the resonance 
domain, i.e. inside a circle of radius T; around Eqi, the elements of the current scattering matrix S vary strongly with 
the energy. In turn, T(E) has also an important variation for the real energies included in the resonance domain. 
Thus the resonances appear usually as peaks in the tunneling coefficient and can be directly seen in the transport 
properties of the structure. While in the case of a ID scattering potential the peaks are light asymmetric maxima^, 
for a 2D scattering potential the peak shapes cover all ranges of the Fano lines, from asymmetric maxima through 
"S-typc" Fano lines up to antircsonances. These profiles have been already seen experimentally for example in the 
conductance of a single-electron transistor—. In the next section we demonstrate that the two-dimensional character 
of the scattering potential and the strong coupling of the quantum system to the contacts allow for the transmission 
profiles which are far from Breit-Wigner lines. 

The representation of the S'-matrix in terms of f2, Eq. (|45[) allows for an efficient numerical procedure to determine 
its poles and the resonances. When the quantum system, for example a quantum dot, is coupled to the contacts it 
becomes open, and the real eigencnergies of the closed problem, Ei , migrate in the lower part of the complex energy 
plane, becoming resonant energies, Eqi = E i — iTi/2, I > 1. On the base of this correspondence we fix an energy E\ 
of the isolated dot and determine the resonance energy E$\ as a solution of Eq. in the complex energy plane. 
The matrix i2 contains contributions from all Wigncr-Eiscnbud functions and energies, i.e. xi an d Eh an d from all 
scattering channels, i.e. all matrix elements of K. Thus the resonance energy Eq\ can strongly differ from E\, and 
only in the case of a very low coupling of the dot to the contacts the eigenenergies of the isolated system, E\, can 
properly approximate the real part of the resonance energy. 

The resonance theory presented above is general and can be applied to a variety of structures with a 2D scattering 
potential, regardless if the geometry of the system is Cartesian or cylindrical. The information about the geometry 
is contained only in the Wigncr-Eiscnbud functions and energies. After solving the eigenvalue problem of the closed 
counterpart of the scattering system one can construct the scattering matrix and analyze it without bearing in mind 
the geometry of the system. 

The expression (|45|) of the 5-matrix shows that all matrix elements of S are singular at the resonance energy. 
That means that all transmission coefficients T nn > between different scattering channels have a similar dependence 
on energy around a resonance, and it is enough to analyze the total transmission which is a sum of them in order to 
characterize the resonance. 



III. MODEL SYSTEMS 

Further we analyze the total tunneling coefficient T(E) for a large energy interval in the case of a quantum dot 
isolated inside of a 2DEG (Cartesian geometry) and in the case of a conical quantum dot in a cylindrical nanowire 
(cylindrical geometry). The transmission peaks are directly connected to the resonances and they have different 
profiles depending on the coupling strength between the quantum system and contacts, but also between resonances. 

A. Quantum dot in two-dimensional electron gas 

We consider here a quite simple dot, a square dot, isolated inside a quantum wire by the constant barriers Vq as 
seen in Fig. [21 The smaller barriers Hi and H2 characterize the coupling between the quantum dot and contacts, 
and the strength of this coupling can be varied individually. Although our model allows for an arbitrary form of the 
potential, we have chosen this square dot in order to compare the scattering functions at the resonant energies with 
the eigenfunctions of an isolated dot (Vq, Hi, H2 — > 00). 

For the numerical calculations we have set d x = d y = 50 nm, and the width of all barrier 20 nm. Thus the region 
where the electrons are localized is about 60 x 60 nm. The barrier which isolates the quantum dot inside the quantum 
wire has been taken as H = 0.4 cV, and the potential energy in the aperture regions Hi = H2 = 0.005 cV. In the 
source and drain contacts the potential energy has been considered as the energy reference, V\ = H = eV. The 
Fermi energy of the electrons has been taken as Ep = 12 meV. For this quantum dot there are four open channels, 
Ni(Ep) = N2(Ep) = 4, and the closed channels do not have a significant contribution to the conductance. 

The quantum dot described above is strongly coupled to the source and drain contacts because the potential energy 
in the aperture regions lies under the Fermi energy. For explaining the transport phenomena through the dot it is 
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v. 




d v — 
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FIG. 2: Potential energy in the 2D quantum wire: constant potential energy in the source and drain contacts, Vi ~ V2, and 
position dependent potential energy in the dot-region. The quantum dot is isolated inside the quantum wire by the barrier 
with the height Vq. The coupling between dot and contacts is set by the potential energy in the aperture regions, 14 1 and Via- 
The electrons inside the dot experience the potential energy Vd- 



necessary to take properly into account the open character of the system and to analyze the transport properties in 
terms of the resonances. The eigenstates which characterize the closed counterpart of the open dot have an infinite 
life time and can not explain the broaden peaks which are experimentally measured in the conductance of an open 
quantum &o&. 

The resonance energies of the considered dot, solutions of Eq. (|49|) . are presented in Fig. [3l For comparison, the 
Wigner-Eisenbud energies, i.e. the eigenenergies of the corresponding isolated dot are also given. Due to the coupling 
of the quantum dot to the contacts the resonance energies migrate in the lower part of the complex energy plane 
and have different widths. There are very narrow resonances associated with the modes of the dot, which are not so 
strongly coupled to the contacts, and broad resonances, which describe modes strongly perturbed by the interaction 
with the reservoirs. For a better understanding of the resonance modes we will examine the localization probability 
distribution density of the electrons for the energies given by the real part of the resonance energy E i . 
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FIG. 3: Resonance energies Eqi (empty circles) of the open quantum dot given in Fig. [2] and the real eigenenergies (filled 
triangles) of the isolated counterpart of the considered quantum dot. 

The potential energy felt by the electrons inside the dot is Vd- This energy can be modified continuously by varying 
the voltage of a plunger gat e 20 ' 22 ' 23 , and the conductance through the dot is measured as a function of Vd- In the 
linear regime experiments, i.e. small source drain biases, and for very low temperatures the conductance G can be 
directly connected to the total tunneling coefficient at the Fermi energy^ i 22 ' 32 , 

G(V d ) = ^-T(E F ;V d ). (50) 

In the above relation the potential energy Vd appears as a parameter in the expression of the total tunneling coefficient. 
A new value of Vd means a new scattering potential and a new scattering matrix. But there is no analytic dependence 
of T on Vd, so that the scattering matrix and after that the conductance have been numerically computed for each 
value of Vd- In Figs. [Hand [5] the conductance is plotted as a function of Ep — Vd- This fact has the great advantage 
that the position of the maxima in conductance are given with respect to the bottom of the quantum well, and in 
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FIG. 4: (Middle part) Conductance (solid line) as a function of the potential energy in the dot region, Vd- The vertical dashed 
lines give the position of the potential energy Voi, I = 1,3,4. (Upper part) The space dependence of the electron probability 
distribution density P n {x,y)/P ma x, P n (x,y) = \^ i n ) (E; x, y)\ 2 , P max = max[P n (i, y)\, n = 1 or n = 2, for E = E F and 
Vd = Voi, I = 1,2,3,4. For the peaks I = 1,2 we have considered the channel number n — 1, while for the peaks I — 3,4, 
n — 2. (Lower part) Resonance energies in the complex energy plane with the real part around Ef- The vertical dashed lines 
correspond to the Fermi energy. 



this way a direct comparison with an infinite quantum well is possible. There are narrow and broad peaks in the 
conductance and in order to understand why they have different profiles we have also plotted the electron probability 
distribution density P n (x,y) = \ijjn (E;x,y)\ 2 for E = Ep and Vd = Voi, f° r the eight peaks considered here, I = 1,8. 
In principle the potential energy Voi is associated with the maximum of the conductance peak, but the conductance 
curve shows also a "S-type" Fano line, and a rigorous method to fix Voi is necessary. 

Using the R-matrix representation of the S-matrix, Eq. (|45[) we can provide an approximative relation for T(E F , Vd) 
around a resonance 

T(E F , V d ) = T(E F , Voi + SV) ~ T(E F - SV, Voi) (51) 

where Voi is the value of Vd for which the real part of the resonance energy (Eoi = Eoi — iTi/2) matches the Fermi 
energy, Eoi = E F , and SV is a small variation with values in the interval (— r;,rj). For a detailed discussion of this 
approach see Ref.— , Appendix A. The expression (f5Tj) of the total tunneling coefficient allows for a direct connection 
to the resonances. We can simultaneously plot the conductance as a function of E F — SV, SV G (— T^T;) and the 
resonances with the real energies in the interval (E F — Ti,E F +Ti). These plots are given in Figs. [4] and [5] The 
dashed vertical lines correspond in each picture to SV = 0, i.e. Vd = Voi in the plots in the middle part and E = Ep 
in the plots in the lower part. From the simultaneous analysis of these graphics it is evident that we can associate 
each peak in the conductance with a resonance I. At the resonance energy the electrons show a strong localization in 
the dot region as seen in the upper part of the figures. A narrow peak corresponds to a resonance energy with a very 
small imaginary part and to a resonance state that is almost decoupled from the contacts, i.e the electron probability 
distribution density is nearly zero in the aperture regions. These are the resonances denoted by (3), (4), (6) and 
(8) in Figs. [5] and [5] The other peaks are broaden and they correspond to resonances with a larger imaginary part 
and to states which are strongly coupled to the contacts. Generally, the peaks for which the resonance states have a 
maximum at y = couple strongly to the contacts and become broaden. These are the peaks (1), (2), (5) and (7) in 
Figs. 2] and E] But there is an exception which does not depend on the parameters of the system: the peak denoted by 
(6) in Fig. [5] Although the probability distribution density has a maximum in the central region of the dot, this state 
is strongly localized inside the dot, and the corresponding peak is very narrow. But this peak is not an asymmetric 
maximum anymore, it has a "S-type" Fano line shape. This behavior can be explained only taking into account the 
interaction between resonances. The state which corresponds to the peak (5) has three maxima in the x-direction 
and one maximum in the y-direction, while the peak (6) has one maximum on x- and three maxima on y-dircction. 
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FIG. 5: (Middle part) Conductance (solid line) as a function of the potential energy in the dot region, Vd- The vertical dashed 
lines give the position of the potential energy Voi, / = 6,8. (Upper part) The space dependence of the electron probability 
distribution density P n (x,y)/ P max , P n {x,y) = \ipn (E; x, y)\ 2 , P max = max[P n (x, y)] for E = E F and V d = V ot , I = 5,6,7,8. 
For the peaks I = 5, 6, 7 we have considered the channel number n = 1, while for the peaks I = 8, n = 2. (Lower part) Resonance 
energies in the complex energy plane with the real part around Ef- The vertical dashed lines correspond to the Fermi energy. 



They are states with the same symmetry in the both directions and they influence each other. As a result, there 
are two hybrid modes, one of them very strongly coupled to the contacts and the other one almost isolated. This 
interaction between resonant states with the same symmetry in the lateral direction is a general phenomenon which 
has at the origin the scattering between different energy channels due to the nonseparable character of the scattering 
potential. For systems with an effective ID scattering potential the interaction between resonances is weak and the 
strong asymmetric Fano line shapes ("S-typc" or antircsonance) do not appear. Also the hybrid modes do not exist 
in this case. 

Besides the potential energy in the dot region Vd, the confinement potential of the dot Vq can be also modified 
by the top gates. Decreasing the strength of the confinement the coupling of the resonance states to the contacts 
increases, and the conductance peaks become broader. We present in Fig. [5] the evolution of the conductance peaks 
(7) and (8) when the potential Vq decreases. The width of the two peaks increases, and the maximum of each peak 
migrates to higher energies. But the shift in energy is different. The state (7) which is strongly coupled to the contacts 
in the aperture regions is much less influenced by the variation of Vq. The state which corresponds to the peak (8), 
with nearly zero probability distribution density in the region of the apertures, can have a coupling to the contacts 
only in the case of a rather transparent confinement barrier Vq. This explains the significant broadening of the line 
shape and also the larger shift of the peak energy. 



B. Conical quantum dot inside a cylindrical nanowire 



We consider a conical quantum dot, embedded in an infinite cylindrical nanowire with the same radius, as is sketched 
in Fig. 7(a) The parameters considered are the height of the dot h = 5nm, the radius of the nanowire R = 5nm 
and the effective mass m* = 0.19mo. We set in our computations d\\ = d z — 16nm and the total number of channels 
(open and closed) N = 8. In our calculations, the results do not change if more channels are added. 

Depending on the band-offsets between the dot material and the host material the potential produced by the dot 
can be repulsive, yielding a quantum barrier, or attractive, yielding a quantum well. We consider here that the dot 
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FIG. 6: Conductance as a function of the potential energy in the dot region, Vd, for different values of the lateral confinement 
potential, Vq. 



yields an attractive potential V(z, r), represented in Fig. |7(b)| by a quantum well of depth Wb = — 0.125eV. 

The total tunneling coefficient versus the incident energy E is plotted in Fig. [8] for different magnetic quantum 
numbers to. The transmission increases with a unity, every time a new channel Ej™^ becomes available for transport, 
i.e. becomes open. The length of the plateaus is given by the difference between two successive transversal mode 
energies, and this length increases with m. Due to the presence of the quantum well, deviations appear from the 
step-like transmission. Just before a new channel gets open there is a dip, i.e. sharp drop, in the tunneling coefficient. 
These dips are owing to modification of the tunneling coefficient due to the evanescent (closed) channels^ 3 .. This is 
a multichannel effect that was also put in evidence in Cartesian coordinates for quantum wires tailored in a two- 
dimensional electron ga a 33 i 34 ' 35 i 36 ' 37 . 

The dips can be understood considering the simple couple-mode mode l 31 ' 33 i 34 i 35 . For a dot surrounded by the host 
material, the scattering potential V(z, r) is not anymore separable, so that the scattering mixes the channel a 31 ' 33 i 34 ' 35 . 




FIG. 7: (a) Sketch of a conical quantum dot embedded into a nanowire with the same radius. The dot yields an attractive 
potential V(z, r), represented in (b) by a quantum well of depth Wb = — 0.125eV. 
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FIG. 8: Total tunneling coefficient as a function of incident energy E for the scattering potential represented in Fig. |7(b)| 
for different magnetic quantum numbers m (continuous line for m = 0, dashed line for |m| = 1) and for the well depth 
W b = -0.125eV. 



As soon as the scattering potential is attractive, the diagonal coupling matrix clement 

V nn (z)= [ (l> n (r)V(z,r)<l> n (r)rdr <0 (52) 
Jo 

acts for every channel n as an effective one-dimensional (ID) attractive potential^!, which always allows for at least 
one bound stat o 46 ' 47 below the threshold of the continuum spectrum. By mixing the channels, this bound state 
becomes a quasi-bound state or resonance, i.e with complex energy, whose real part gets embedded into continuum 
spectrum of the lower channel and the imaginary part describes the width of the resonance. These resonances can 
be seen now as dips in the tunneling coefficient. The energy difference between the position of the dips and the 
next subband minima Ej"£ gives the quasi-bound state energy. The positions of the dips, i.e. the quasi-bound state 
energy, depend on the channel number n and on the magnetic quantum number m and, of course, on the detailed 
system parameters. In Cartesian coordinates the specific symmetry of the channels (odd and even) do not allow for 
dips in the first plateau^. In the cylindrical geometry this symmetry is broken, so that we obtain a dip in front of 
every plateau. Our numerical method allows for a high energy resolution in computing the tunneling coefficient, so 
that we were able to find the dips also in front of the higher-order plateaus. 

Further insight about the quasi-bound states of the evanescent channels can be gained looking at the wave functions, 
whose square absolute value \tpn (E; z,r)\ 2 gives the localization probability density. The R-matrix formalism allows 



us to produce high resolution maps of the wave functions inside the scattering region, see Eq. (|33|) . In Figs. 9(a) 
|9(b)| the localization probability distribution density is represented in arbitrary units, for an electron incident from 
source (s = 1) and with a total energy corresponding to the dips in Fig. [5J The total energy E and the channel 



on which the electron is incident, are specified at every plot. Let discuss Fig. 9(a) The total energy E — 0.199eT^ is 



less than the energy of the second transversal mode, E^ 2 = 0.244eT^, so that only the first channel is open. Thus the 



incident wave from the source contact is node-less in r-direction. But, as it can be seen in Fig. 9(a) the scattering 
wave function inside the scattering region has a node in the r-direction, i.e. position in r where the wave function is 
zero. This means that the wave function corresponds to the quasi-bound state splitting off from the second transversal 
mode, which is an evanescent one. The quasi-bound state is reachable now in a scattering formulation due to channel 
mixing. The scattering wave function has a pronounced peak around the scattering potential, i.e. z 6 [0, 5]nm, and 
decreases exponentially to the left and to the right. On the left side of the scattering potential one observes the 
interference pattern produced by the incident and the reflected waves, while on the right side there exists only the 
transmitted wave. 



The scattering wave function considered in Fig. 9(b) has the energy less than the third transversal channel, 



-Ej° 3 = 0.6006eV, so that the incident part of the scattering state on the second mode n = 2 has one node in re- 
direction. But the scattering function shows inside the scattering region two nodes in the r-direction, so it corresponds 
to a quasi-bound state splitting off from the above evanescent channel, the third one. 

One gets similar pictures for all m-valucs, with the difference that for m ^ the wave functions are zero for r = 0. 
In Fig. [9(a)] and [9(b)1 one can observe that the transmitted part of the scattering wave function is zero, in agreement 
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FIG. 9: Localization probability distribution density, \ipn\E; z,r)\ 2 , for an electron with m — 0, incident from the reservoir 
s = 1 into the channel n and with the total energy E. The energies are the dips in Fig. [8] 



with the resonant backscattcring specific to the quasi-bound states of the evanescent channe b 33 i 34 . Increasing the 
strength of the attractive potential one can see more dip a 31 ' 36 in the tunneling coefficient. Another systems embedded 
inside the cylindrical nanowire, like a cylindrical dot, a quantum ring or a double barrier heterostructure, which also 
show a similar behavior, were studied in Refi 3 -^. 

We have analyzed until now only the classical allowed energy domain, with a continuous, double degenerated 
spectrum. The wave functions of the electrons for these energies are extended states, presented here as scattering 
states. The classically forbidden spectrum contains the bound states or the localized states. The R-matrix formalism 
presented here can provide also these states, as long as the boundary points ±d z are far enough from the quantum 
dot, so that the bound states fulfill the Neumann boundary condition l[25|). In such a way, the energies of the bound 
states are the negative Wigner-Eiscnbud energies and the wave functions for the bound state are the corresponding 
Wigner-Eisenbud functions. For the conical dot presented here, there is only a bound state represented by the lowest 



Wigner-Eisenbud energy, Eb = E\ < and ipb(Eb;z,r) = %i(z,r). In Fig. 10(a) we present the energy value Eb of 
the bound state together with a cut of the potential energy along r = 0. The corresponding wave function, namely 
the absolute value square, is represented in Fig. |10(b)| One can see that the bound state is mainly localized inside 
the quantum dot. 




FIG. 10: The bound state for a conical quantum dot inside the nanowire, as in Fig. [7] (a) Value of the bound energy 
Eb = — 0.013eV (dashed line). The potential energy along r — is represented by the solid line, (b) The absolute value square 
of the wave function corresponding to the bound state, \ipb(Eb', z, r )| 2 . 



IV. SUMMARY AND DISCUSSION 



We have presented a general theory for computing the scattering matrix and the scattering wave functions for a 
general finite-range extended scattering potential in two dimensions. The theory is based on the R-matrix formalism, 
which allows a semi-analytical treatment of the scattering problem, yielding in such a way a powerful and efficient 
numerical method. 
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This formalism was applied to a quantum dot defined inside a two-dimensional electron gas, as well to a conical 
quantum dot embedded inside a cylindrical nanowire. 

It is pointed out the role of the evanescent channels, which for a nonseparable attractive scattering potential in a 
multi-channel nanowire produces resonant dips in the tunneling coefficient. Furthermore, the cylindrical symmetry 
does not yield the same "selection rules" for tunneling coefficient as the Cartesian symmetry. 

It is also presented a general resonance theory, which shows that the two-dimensional character of the scattering 
potential and the strong coupling of the quantum system to the contacts allow for the transmission profiles which 
ranges from asymmetric Fano line shapes, through "S-type" Fano lines until antircsonanccs. 

Detailed maps of localization probability distribution density sustain the physical interpretation of the resonances 
(dips and peaks) found in the studied heterostructures. 



Acknowledgments 

It is a pleasure for us to acknowledge the fruitful discussions with Klaus Gartner, Vidar Gudmundsson, Andrei 
Manolescu and Gheorghe Nenciu. One of us (P.N.R.) also acknowledge partial support from German Research 
Foundation through SFB 787 and from the Romanian Ministry of Education and Research through the Program 
PNCDI2, Contract number 515/2008. 



racec@wias-berlin.de 

roxana@physik.tu-cottbus.de 

neidhard@wias-berlin.de 



* Electronic address 
' Electronic address 
Electronic address: 

Semiconductor Industry Association, International Technology Roadmap for Semiconductors, http://www.itrs.net. 
B. Doyle, R. Arghavani, D. Barlage, S. Datta, M. Doczy, J. Kavalieros, A. Murty, and R. Chau, Intel Technology Journal 
6, 42 (2002). 

T. Bryllert, L.-E. Wernersson, T. Loewgren, and L. Samuelson, Nanotechnology 17, S227 (2006). 
K. H. Yeo et al, Tech. Dig. - Int. Electron Devices Meet. p. 539 (2006). 

K. H. Cho, K. H. Yeo, Y. Y. Yeoh, S. D. Suk, M. Li, J. M. Lee, M.-S. Kim, D.-W. Kim, D. Park, B. H. Hong, et al, Appl. 
Phys. Lett. 92, 052102 (pages 3) (2008). 

A. D. Wieck and K. Ploog, Appl. Phys. Lett. 56, 928 (1990). 

J. Gores, D. Goldhaber-Gordon, S. Heemeyer, M. A. Kastner, H. Shtrikman, D. Mahalu, and U. Meirav, Phys. Rev. B 62, 
2188 (2000). 

H.-S. Wong, K. Chan, and Y. Taur, IEDM Tech. Dig. pp. 427-430 (1997). 

X. Huang, W.-C. Lee, C. Kuo, D. Hisamoto, L. Chang, J. Kedzierski, E. Anderson, H. Takeuchi, Y.-K. Choi, K. Asano, 
et al, IEDM Tech. Dig. pp. 67-70 (1999). 

B. Doyle, S. Datta, M. Doczy, S. Hareland, B. Jin, J. Kavalieros, T. Linton, A. Murthy, R. Rios, and R. Chau, IEEE Electron 
Dev. Lett. 24, 263 (2003). 

J. Xiang, W. Lu, Y. Hu, Y. Hu, H. Yan, and C. M. Lieber, Nature 441, 489 (2006). 

M.T.Bjork, B.J.Ohlsson, C. Thelander, A. Persson, K. Deppert, L. Wallenberg, and L. Samuelson, Appl. Phys. Lett. 81, 
4458 (2002). 

J. Wensorra, K. M. Indlekofer, M. I. Lepsa, A. Forster, and H. Liith, Nano Lett. 5, 2470 (2005). 

F. Qian, Y. Li, S. G. Caronak, H.-G. Park, Y. Dong, Y. Ding, Z. L. Wang, and C. M. Lieber, Nature Mater. 7, 701 (2008). 
Y. Hu, H. O. H. Churchill, D. J. Reilly, J. Xiang, C. M. Lieber, and C. M. Marcus, Nature Nanotechnol. 2, 622 (2007). 
B. Tian, X. Zheng, T. J. Kempa, Y. Fang, N. Yu, G. Yu, J. Huang, and C. M. Lieber, Nature 449, 885 (2007). 
K. H. Cho et al., Tech. Dig. - Int. Electron Devices Meet. p. 543 (2006). 

D. Vasileska, D. Mamaluy, H. R. Khan, K. Ravela, and S. M. Goodnick, J. Comput. Theor. Nanosci. 5, 999 (2008). 
L. Smrcka, Superlatt. and Microstruct. 8, 221 (1990). 

U. Wulf, J. Kucera, P. N. Racec, and E. Sigmund, Phys. Rev. B 58, 16209 (1998). 

E. Onac, J. Kucera, and U. Wulf, Phys. Rev. B 63, 85319 (2001). 
E. R. Racec and U. Wulf, Phys. Rev. B 64, 115318 (2001). 
P. N. Racec, E. R. Racec, and U. Wulf, Phys. Rev. B 65, 193314 (2002). 

G. A. Nemnes, U. Wulf, and P. N. Racec, J. Appl. Phys. 96, 596 (2004). 
G. A. Nemnes, U. Wulf, and P. N. Racec, J. Appl. Phys. 98, 084308 (2005). 
T. Jayasekera, M. A. Morrison, and K. Mullen, Phys. Rev. B 74, 235308 (pages 9) (2006). 
U. Wulf, P. N. Racec, and E. R. Racec, Phys. Rev. B 75, 075320 (2007). 

28 J. Behrndt, H. Neidhardt, E. R. Racec, P. N. Racec, and U. Wulf, J. Differ. Equ. 244, 2545 (2008). 

29 G. Mil'nikov, N. Mori, and Y. Kamakura, Phys. Rev. B 79, 235337 (pages 5) (2009). 

30 E. P. Wigner and L. Eisenbud, Phys. Rev. 72, 29 (1947). 

31 P. N. Racec, E. R. Racec, and H. Neidhardt, Phys. Rev. B 79, 155305 (2009). 



18 



E. R. Racec, P. N. Racec, and U. Wulf, unpublished. 
R F. Bagwell, Phys. Rev. B 41, 10354 (1990). 

S. A. Gurvitz and Y. B. Levinson, Phys. Rev. B 47, 10578 (1993). 
J. U. Nockel and A. D. Stone, Phys. Rev. B 50, 17415 (1994). 

J. H. Bardarson, I. Magnusdottir, G. Gudmundsdottir, C.-S. Tang, A. Manolescu, and V. Gudmundsson, Phys. Rev. B 70, 
245308 (2004). 

V. Gudmundsson, Y.-Y. Lin, C.-S. Tang, V. Moldoveanu, J. H. Bardarson, and A. Manolescu, Phys. Rev. B 71, 235302 
(2005). 

B. G., B. J. A., and F. R., Sol. State Phys. 44, 229 (1991). 

F. F. Fang and W. E. Howard, Phys. Rev. Lett. 16, 797 (1966). 

M. Baro, H.-C. Kaiser., H. Neidhardt, and J. Rehberg, Rev. Math. Phys. 16, 281 (2004). 

H. Schanz and U. Smilansky, Chaos Solitons & Fractals 5, 1289 (1995). 

M. Luisier, A. Schenk, and W. Fichtner, J. Appl. Phys. 100, 043713 (2006). 

A. M. Lane and R. G. Thomas, Rev. Mod. Phys. 30, 257 (1958). 

M. Bxittiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. Rev. B 31, 6207 (1985). 

A. Bohm, Quantum Mechanics (Springer, New York, 1993). 

B. Simon, Ann. Physics 97, 279 (1976). 
M. Klaus, Ann. Physics 108, 288 (1977). 



